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Abstract 

We build upon work by C. K. Law [Phys. Rev. A 71, 034306 (2005)] to show in general that 
the entanglement between two fermions largely determines the extent to which the pair behaves 
like an elementary boson. Specihcally, we derive upper and lower bounds on a quantity Xn+i/xn 
that governs the bosonic character of a pair of fermions when N such pairs approximately share 
the same wavefunction. Our bounds depend on the purity of the single-particle density matrix, 
an indicator of entanglement, and demonstrate that if the entanglement is sufficiently strong, the 
quantity Xn+i/xn approaches its ideal bosonic value. 
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Under what circumstances can a pair of fermions be treated as an elementary boson? 
Many authors have done detailed studies of this question, as it applies, for example, to 
atomic Bose-Einstein condensates , excitons and Cooper pairs in superconductors 
jsj. In a 2005 paper, C. K. Law presented evidence that the question can be answered in 
general in terms of entanglement: two fermions can be treated as an elementary boson 
if they are sufficiently entangled jsl. Consider, for example, a single hydrogen atom in 
a harmonic trap. Within the atom, the proton and electron are strongly entangled with 
respect to their position variables; for example, wherever the proton might be found — it 
could be anywhere in the trap — the electron is sure to be nearby. Law suggests that this 
entanglement is the essential property underlying the (approximate) bosonic behavior of the 
composite particle, allowirig, for example, a collection of many hydrogen atoms to form a 
Bose-Einstein condensate |7|. 

Specifically, his hypothesis can be expressed as follows: For a single composite particle 
in a pure state, let P be the purity of the reduced state of either of the two component 
fermions — P is small when the entanglement between the two particles is large (see below 
for the definition) — and let be the number of composite particles that approximately 
share the given quantum state. Then the composite particles can be treated as bosons as 
long as NP ^ L That is, according to this hypothesis, the quantity 1/P roughly quantifies 
the number of particles one can put into the same pure state, before the composite nature 
of the particles begins to interfere appreciably with their ideal bosonic behavior. 

Law's argument in support of this hypothesis assumes a two-particle wavefunction within 
a certain class, characterized by a specific form of the eigenvalues of the reduced density 
matrix of either particle, and he notes that it would be desirable to extend the argument to 
more general wavef unctions. Such a generalization is the aim of the present paper. With no 
restrictions on the form of the two-particle wavefunction, we use the purity to place upper 
and lower bounds on Law's measure of bosonic character, and we show that these bounds 
are the tightest possible of the given form. In this way we obtain a more general connection 
between entanglement and bosonic character. 

Before proceeding to our general argument, it may be instructive to consider the special 
case of the hydrogen atom. Let \E'(-R, r) be the wavefunction of a single hydrogen atom in 
a harmonic trap, with R and r being the position coordinates of the proton and electron, 
respectively. For simplicity we assume that the proton is sufficiently massive compared to 
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the electron that we can write this wavefunction as 

<i/{R,r) =ilj{R)(f){f- R), (1) 
where ip is the ground-state harmonic oscillator wavefunction 

^(R) = :;^2-M-RV2b% (2) 

and is the ground-state wavefunction of the electron in a hydrogen atom: 

exp(-r/ao). (3) 



Here Cq is the Bohr radius and 6 is a length parameter characterizing the size of the trap. 
The purity P of the reduced state of either of the two particles is defined by 

P = Trp2, (4) 

where p is the density matrix of the particle in question. (Because the pair is in a pure 
state, the purities of the two particles are guaranteed to be equal.) Note that P takes values 
between zero and one. For the hydrogen atom, the purity of the proton is given by 

P = j p{R,R')p{R',R)dRdR', (5) 

where the proton's density matrix p is 

p{R,R') = j ^{R,r)^*{R\f)df. (6) 

Inserting Eq. ([T]) into the definition of P, we find that 

P = J\^{R)\^\ij{R')\^\a{R' -R)\^dRdR' (7) 

= jmRrmR + q)mq)?dRdq, (8) 

where 

'^(q) = / 4'{r)(t)*{r — q) dr. (9) 



But the range of a is comparable to the Bohr radius and much smaller than the dimension 
of the trap. So we can reasonably replace \ip{R + in Eq. (jH]) with and write 

P = J \^{R)\UR j \a{q)\^dq. (10) 



The integrals can be done, and one finds that 

P=J^f-_^)\ (11) 

Thus the purity depends, not surprisingly, on the ratio of the volume of an atom to the 
volume of the trap, and Law's condition NP <^ 1 essentially says that the space available 
to each atom must be large compared to its size. This condition is in rough agreement with 
the condition that the number of atoms be small compared to the maximum occupation 
number as computed in Ref. 

We now turn to the general argument. 

Consider a composite particle formed from two distinguishable, fundamental fermions 
A and B with wavefunction \E'(xyi,XB). (Here the x's could be vectors in any number of 
dimensions.) Writing this wavefunction in its Schmidt decomposition yields 

^{xa,xb) = $^Af 0i^)(x^)0f )(xb). (12) 
p 

Here (pp^^ and (pp^^ are the Schmidt modes, constituting orthonormal bases for the states 
of particles A and B, and the Ap's, which are the eigenvalues of each of the single-particle 
density matrices, are nonnegative real numbers satisfying Ylp^p ~ ^- ■'^^ terms of the Ap's, 
the purity can be written as 

p 

Again, a small value of the purity indicates a large entanglement. 

In terms of creation operators, the state \1/(x^,xb) can be written as 

^f{xA,XB) = J2K^'<K\^)^ (14) 

P 



where aj^ creates an A particle in the state (f)p (xa), bj^ creates a B particle in the state 
(Pp^Hxb), and |0) is the vacuum state. The composite particle creation operator c^, which 
creates a pair of A and B particles in the state \l/(a;^, x^), is defined to be 

^' = E^r* (15) 
p 

Our analysis, like Law's, aims to determine to what extent the operators and c act 
like bosonic creation and annihilation operators when applied to a state consisting of 
composite particles. 



Consider the state obtained by antisymmetrizing the product state 
a;^"*) ■ ■ ■ \[^(a;^'', x^^). In terms of the creation operator c^, we can write the 
properly antisymmetrized state as 



\N) = 4^x-^'%r\o). 



(16) 



Here is a normahzation constant necessary because is not a perfect bosonic creation 

nn 

operator. The quantity Xn given by |6|, 18| 



X7v = ^(0|c^(ct)^|0)= Yl VAp....Ap,. 



(17) 



pi ■ • ■ PJV 
all different 



(This expression gives xn = if the number exceeds the number of Schmidt modes with 
nonzero Schmidt coefficient. In that case (c^)^|0) =0 and we cannot define the state |A^).) 
For an ideal boson, we would have xn = 1- 

Note that c^N) is not necessarily equal to \/ N + I + 1). Rather, it follows from the 
definition f lT6|) that 

c^\N) = aN+iVNTl\N + l), (18) 

where 

^ (19) 



On 



Xn 



Xn-1 



Similarly, instead of c\N) = vNlN — 1), we have 



c\N) = aWN\N -l) + \eN), 



(20) 



where |eAr) is orthogonal to — 1). For perfect bosons, we would have {cnI^n) = 0, but 
the actual value is p, 



(eAr|eAr) 



^_XN+1_^ 



Xn 



Xn+1 



(21) 



Xn \Xn~i Xn 

In the Appendix, we show that the ratio Xn+i/Xn which appears in Eqs. (fT9|) and ( 12T|) 
is strictly non-increasing as increases (more precisely, we show that Xn ~ Xn+iXn~i is 
non-negative), so that the quantity in parentheses in Eq. (12 ip is non-negative. It follows 
that both otv and (eisfl^N) will be within a small amount 6 of their bosonic values when 
Xn+i/xn > I — S. One can also show jsl, Q] that 



Xn J 



(22) 



which is within 2S of its ideal bosonic value, 1, under the same condition. We therefore 
follow Law in using the ratio xn+i/xn — we call it the "xAf-ratio" — as our indicator of 



bosonic character \6 



One might wonder why we confine our attention to quantities involving only the state 
|A^) and nearby states, rather than insisting that the operator c act like a bosonic operator 
on the whole subspace spanned by {|0), . . . , |A^)}. The reason is that we are interested in 
a state that approximates |A^), and we wish to quantify the degree to which the system 
behaves like a collection of bosons when a composite particle is added to or removed from 
this state. Hence our focus on Xn+i/xn as the quantifier of bosonic character rather than 
Xn itself. We note that because the XAf-ratio is non-increasing with A^, a lower bound on 
XAT+i/X-M-Will also be a lower bound on xn'+i/xn' for all A^' < A^. However, as one can see 
in Ref. |8|, this fact is not sufficient to guarantee that xn itself is close to unity whenever 
Xn+iIxn is. 

In the remainder of the paper we prove two inequalities relating the Xiv-ratio to the 
purity. 

The first is a lower bound: Xn+i/xn > 1 — NP. To show this, we consider the quantity 
Xn+i — Xn{X ~ NP) and show that it must be non-negative. 



XN+i -xn{i-np) 



ApiAp2 . . . Ap^_^^ ~ I 1 ~ ^ X] -^p 1 Ap,Ap2 . . . Apj^ 

...pjv+i \ p J Pi-.-Pn 



Pi ■ 

all different all different 



2 



~ ApiAp2 . . . Ap^_^^ - ^ Xp^Xp^ . . . Xp^Xp^_^^ + N ^ ApiAp2 . . . Ap^A 

pi...p]V+l pi.-.pN pi.-.pN 

all different all different; all different; 

pjv+i free Piv+i free 

(23) 

Note that the first two sums of the last line have many terms in common, which therefore 
cancel out. The only terms remaining from those sums are the terms in the second sum for 
which the value of pn+i is equal to the value of one of the indices pk with k = 1, . . . ,N. 
Each of these N possibilities yields the same result; so we can combine those first two sums 
into the expression 

E ^l^p^---^.- (24) 

Pi • ■ ■ Pjv 
all different 
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We therefore have 



XN+i-XN{i- NP) - -N ^ Xl^Xp2- ■ ■ Xpj^Xp^^^ + N ^ W • • • -^Piv^Piv+r 



P1...PJV P1---PN 
all different; all different; 

pjv+i free PJV+i free 



(25) 

Again the two sums have many terms in common. Cancelhng these terms leaves 



Xn+i-Xn{1-NP) = N{N-1) XlXp2...Xp,-N{N-l) Yl 

pi-.-PN p\...pN 
all different all different 

= N{N - 1) Y ^P^^P^ ■■■^PM - W) • 

P1---PN 

all different 

(26) 

Now, Eq. fl^Ul) can be rewritten as 



Xn+1 - Xn{'^ - NP) - — ^— ^ Y ■ ■ ■ Kn {^li + K2 ~ '^^Pi^Pi) 

pi-.-PN 

all different 
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Pi ■■■PN 
all different 

thus yielding the bound 

>1-NP. (28) 

Xn 

This bound shows that a sufficiently small purity entails nearly bosonic character as 
quantified by Xn+i/xn- We now derive a bound in the other direction, showing that a 
nearly bosonic value of Xn+i/xn implies a small purity. For this purpose we start with 

(1 - P)Xn - Xn+1 = I 1 ~ I Xpj^Xp2 . . . Xp^ - ^ ApiAp2...A 



Vp2 . . . ^vpjv^i 

all different all different 



y^ ^Pl^P2 ■ ■ ■ ^PN^PN+1 y^ -^Pl-^P2 • • • -^PJV-^PjV+1 y^ -^Pl-^P2----^PAr+l 

P1...PJV pi... PAT pi...piV+l 

all different; all different; all different 

Pjv+1 free pjv+i free 

(29) 
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By combining sums as before (inserting the identity 1 = X]p^r+i '^pn+i when needed), we get 



P1---PN+1 P1---PN 
all different all different 



(30) 

Combining this result with our earher inequahty (Eq. (j2H|) ). we have 

l-NP< < 1 - P. (31) 

Xn 

We have thus put upper and lower bounds on the XAf-ratio of a composite particle made 
of two distinguishable fermions, in terms of the entanglement of the pair. We have not 
specified anything about the form of the wavefunction of the composite particle; so the link 
between the XAf-ratio and entanglement is established in general. 

The lower bound in Eq. (13T|) is in fact as strong a bound as one could hope to derive 
in terms of purity, in that the bound is achievable: if there are M nonzero Schmidt modes 
and \p = 1/M, then, by Eq. (^Tij, Xn+iIxn = 1 — NP as long as is less than M. This 
lower bound is also achieved by wavefunctions in the class Law considers — this class includes 
double Gaussian wavefunctions — in the limit NP ^ 1. Because Eq. (l30ll is never zero unless 

= 1 (in which case it is always zero), our upper bound is not, for general A^, achievable. 
Nevertheless, it is the best possible upper bound of the form Xn+i/xn < 1 — bP, whether 
or not h depends on A^. This is because for any value of h greater than 1, there exists a 
distribution of Schmidt coefficients that makes 1 — hP negative — it suffices to make one of 
the coefficients very large — whereas Xn+i/xn is certainly non-negative. We note also 
that there can be no upper bound of the form 1 — hP^ with r less than 1, because such a 
bound would contradict our lower bound when P is small. 

We have considered in this paper only a single wavefunction '^{xaiXb) of the composite 
particle. One would also like to investigate whether, for several orthogonal wavefunctions 
\E'j(x^,XB), the corresponding creation operators c] approximately satisfy the bosonic rela- 
tion [cj, c^] = for j 7^ k. (The relation [cj, c^] = will automatically be satisfied because of 
the anticommutation of the underlying fermionic operators.) If the relevant deviation from 
this commutation relation similarly diminishes to zero as the entanglement of each wave- 
function increases, one will then have further evidence for the proposition that entanglement 
is crucial for determining whether a pair of fermions can be treated as a boson. 
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Taking this idea to its logical conclusion, Law notes that two particles can be highly 
entangled even if they are far apart. Could we treat such a pair of fermions as a composite 
boson? The above analysis suggests that we can do so. However, we would have to regard 
the pair as a very fragile boson in the absence of an interaction that would preserve the pair's 
entanglement in the face of external disturbances. On this view, the role of interaction in 
creating a composite boson is not fundamentally to keep the two particles close to each 
other, but to keep them entangled. 

We thank Fred Strauch, as well as the participants of the quantum information group at 
Perimeter Institute, for their stimulating questions and valuable suggestions. Research at 
Perimeter Institute is supported by the Government of Canada through Industry Canada 
and by the Province of Ontario through the Ministry of Research and Innovation. 

APPENDIX: PROOF THAT xlf - Xn+iXn-i > 

Let us use the symbol Yl' to indicate a sum over all the indices appearing in the summand, 
with the restriction that they must all take distinct values. We have, then, 

Xn ~ XN+iXN-1 — Yl'^ri ■ • ■ Kn Y'Ki ' ' ' ~ Yl'^ri ' ' ' Kn+i Y! Ki ' ' ' Kn-i ■ (32) 

We now treat separately the sum over Pat in the first term and the sum over r^+i in the 
second term, obtaining 

Xn~ XN+iXN-1 = Z^'-^ri ■ ■ ■ Kn [ ' ' ' Kn-1 ~ Y'Ki ' ' ' Kn-1 {Ki + • • • + '^pjv-i)] 

~ [Y ^ri ' ' ' Kn ~ Y Ki' ' ' Kn {Ki + • ' • + Kn)] Y Ki ' ' ' Kn-i (33) 
— Y Ki'-- Kn Y Ki'" Kn-1 [Ki H \- Kn ~ {Ki ^ '^Pjv-i)] ■ 

We now separate the sum over the p's into two parts: (i) the part in which pi, . . . ,Pn-i have 
the same restrictions as ri, . . . , tat-i, and (ii) the part in which one of the p's has the same 
value as Tat. So the expression becomes 

YKn Z^^rjv'^^i ■ ■ ■ '^'■Jv-i Y^rN'^Pi ' ' ' Kn-1 [Kn + {Ki H h Kn-i) ~ {Ki ^ Kn-i)] 

+ {N - 1)Y'^In Y^^TN^'-i ' ' ' '^rN-l Y'^TN^P^ ' ' ' \n-2 [{Ki H ^ Kn-i) " (V ^ -^PiV-2)] • 

In the first line, everything in the square bracket vanishes except for Xj.^, because of the 
symmetry between the r's and the p's. Thus the first line is non-negative. Meanwhile the 
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expression in the second line that appears within the sum over r^v has the same form as the 
last line of Eq. flHHj) . except with one fewer r index and one fewer p index. Therefore the 
same maneuver can be repeated again and again, until at last we are left with an expression 
that is manifestly non-negative, because there are no more of the negative p terms. It follows 
that 

Xn - Xn+iXn-1 > 0. (34) 
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